The Non-Equilibrium Reliability of Quantum Memories 



0^ 
O 
O 
(N 

X) 
m 

a.; 

a- 
^; 



en 
> 

00 

p 

O 
oo 

o 



X 



Alastair Kay^'^'^ 

^ Max-Planck-Institut fiir Quantenoptik, Hans-Kopfermann-Str. 1, D- 85748 Garching, Germany 

^Gentre for Quantum Gomputation, DAMTP, Centre for Mathematical Sciences, 

University of Gambridge, Wilberforce Road, Cambridge GB3 OWA, UK 

''Centre for Quantum Technologies, National University of Singapore, 3 Science Drive 2, Singapore 117543 

(Dated: February 23, 2009) 

The ability to store quantum information without recourse to constant feedback processes would 
yield a significant advantage for future implementations of quantum information processing. In 
this paper, limitations of the prototypical model, the Toric code in two dimensions, are elucidated 
along with a sufficient condition for overcoming these limitations. Specifically, the interplay between 
Hamiltonian perturbations and dynamically occurring noise is considered as a system in its ground 
state is brought into contact with a thermal reservoir. This proves that when utilizing the Toric 
code on A^^ qubits in a 2D lattice as a quantum memory, the information cannot be stored for a 
time 0{N). In contrast, the 2D Ising model protects classical information against the described 
noise model for exponentially long times. The results also have implications for the robustness of 
braiding operations in topological quantum computation. 



Introduction: An essential element to the rapidly 
evolving technologies involved in the implementation of 
quantum information processing will be the ability to 
store quantum states in the presence of noise. The notion 
of reliable storage of a qubit requires that there should 
be a family of storage solutions parameterised by their 
size, N, and that the time for which the information sur- 
vives should grow exponentially in N. That such solu- 
tions exist is guaranteed by the theory of fault tolerance 
[l|, which states that provided errors occur locally, and 
only at an error rate below a threshold rate which is in- 
dependent of N, suitable encodings of the data can be 
found, although regular error correction cycles are re- 
quired, with massive overheads. 

Evidently, a device more akin to a quantum hard drive 
is desirable, where the data can be input and later read- 
out without intervention in the meantime. Without any 
error correction whatsoever, this idea is clearly flawed; 
any state |^) of iV qubits can be converted into an or- 
thogonal state by a local rotation on any single qubit. To 
see this, consider the Schmidt decomposition of |^) on 
a particular qubit i, |\I>) = a |0)j |V'o) + Z? |l)i I'fpi), where 
|a|2 -t- |/3|2 = 1 and (V'olV'i) = 0. A bit-flip, X, on this 
qubit has the effect (\l/|Xi |^) = 0. For a noise model 
with a constant per-qubit error rate, the length of time 
for which j^P) survives decreases with A^. 

The intermediate regime involves encoding a qubit 
state within a subspace (or subsystem) of N qubits, and 
allowing a single round of error correction when that state 
is read-out. Reliable storage remains feasible since the 
logical gate operations Xl and Zl which act on these 
states, \0l) and |1l), can be non-local, so that local ro- 
tations do not mix the subspace, and can thus be detected 
and corrected at read-out. One concept to ensure that 
these local errors do not build up is to create a Hamil- 
tonian in which the states \0l) and \1l) are degenerate 
ground states. If this Hamiltonian has an energy gap A 
up to the next excited state, then one might hope that 
for any temperature ksT <C A, thermal noise would be 



unable to mix the ground states. This motivated Kitaev 
'2\ to introduce the Toric code in two dimensions. 

Following the introduction of a Hamiltonian H, one 
must consider both noise due to interaction with an en- 
vironment and also perturbations in H. In this paper, 
we consider the interaction between these two types of 
error to prove that the Toric code in 2D does not provide 
reliable storage of information, and include the necessary 
details for the generalisation of the technique to other sta- 
bilizer Hamiltonians. These reveal a sufficient condition 
(for both stabilizer and non-stabilizer Hamiltonians) such 
that the presented test is passed, which requires models 
to exhibit a property known as string tension. The 2D 
Ising model displays this for classical data, as does the 
Toric code in 4D for quantum data [3|. 

The Toric Code in two dimensions is the paradigmatic 
construction of a quantum memory. Consider a plane 
with orthogonal vectors x and z and periodic boundary 
conditions. Qubits are located at positions 2ix + 2jz and 
(2i + l)x + (2j + l)z for all integers < i, j < A^ - 1, and 
they interact via a Hamiltonian 
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Zi,j + Xi . 



X2i,2jX2i+l^2j + lX2i 



2i+l,2j+lZ2i+2,2jZ2i+l,2j-l, Xij — 

is the stan- 



i,2j+iX2i,2j+2 and Aij 
dard Pauli A-operator acting on the qubit at position 
ix + jz. The ground state is 4-fold degenerate, encoding 
two logical qubits, and the loops Z^i ~ Yli ^2j+i,2j+i 
and Zl^2 — Yii ^2i,2j define the Z-basis of the two qubits 
they encode. In order to convert between these states, a 

string of X operators such as TTi^o -^^i,^] needs to be 
applied. It is typically argued 2] that any perturbation 
6H = S^^hi, where hi acts on a constant number of 
qubits in a local way, and satisfies \\hi\\ < 1, only splits 
the degeneracy by an amount S , and thus information 
can survive for a time (5~^. This is because a product 
of 0{N) of the hi is required to describe a logical rota- 



tion within the degenerate subspace, and hence 0{NY^ 
order perturbation theory is required [Ij|. Evidence is 
mounting [1, H, S [1, 0] that the Toric code in 2D is not 
protected in the same, exponential, way against interac- 
tions with a thermal environment, a property known as 
self-correction. It is this which we wish to prove. The 
noise models that have so far been explicitly discussed 
have been based on a Monte-Carlo approach 3] which has 
the advantages that it is easily analyzed, and is known to 
have the thermal state as its steady state. However, there 
is no known physical process that results in this type of 
noise. Alternative treatments mostly consider whether 
storage is possible when the system is in equilibrium with 
a thermal bath [3, 3, [a, B 0| j but are beginning to con- 
sider the time required to approach equilibrium from an 
initial state ^] . All these approaches take it as given that 
Hamiltonian perturbations are not destructive. 

In order to discuss how Hamiltonian perturbations in- 
teract with errors that enter the system, we impose the 
same demands as required for the derivation of fault- 
tolerance - robustness of storage must be with respect 
to an adversarial model i.e. all possible error combina- 
tions subject to some physically motivated restrictions. 
By selecting one particular example, an upper bound on 
the protection given in the worst-case can be derived. 
We avoid a full treatment of an arbitrary environment 
by considering a minimal degree of interaction of such 
an environment as introducing a single fault Xq.o- With- 
out loss of generality, the system is considered to be ini- 
tialised in a ground state jf/') that is the -1-1 eigenstate of 
both Zl^i and Z^ 2- The aim is to use the Hamiltonian 
perturbations to propagate the error through a sequence 
of intermediate operations Ui — Y[i=Q ^2i,o until Un-2 
has been implemented, flipping the logical qubit, up to 
a single local error [i3|. If a Hamiltonian achieves this 
transport in a time i, then data cannot be stored reli- 
ably for longer times. We proceed as in l§\ by defining 
a subspace of interest, spanned by the basis {Ui \4>)} for 
i — 1 ... N — 2. The Hamiltonian for this subspace is de- 
scribed by an (iV— 2) x (iV— 2) diagonal matrix with equal 
diagonal elements (equal to the energy of a single exci- 
tation). There are two perturbations that we can add 
that preserve this subspace structure. The first allows 
alteration of individual diagonal terms, while the second 
is a propagation operator, creating off-diagonal matrix 
elements. 
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i = l 
otherwise 



These terms combine to give a perturbation 

JV-3 

5H = -^5 2_^ >/j(A'2i+i,o(l — ZifiZiJ^ifl)) 



i=0 

Af-2 



\5Y,B^{'^-Z^fl)■ 



The subspace oi H -\- 5H containing the evolution can 
be written as an effective Hamiltonian iJeffj which is 
a tridiagonal matrix with matrix elements {i\Hcs\i) = 
(tM UJEUi IV') = 2A-h5B, and {i\ H^s N + 1) = SJ^. Ma- 
trices of this form have been heavily studied for their 
ability to perfectly evolve a state \i) into |A'' — 2 — i). For 
example, selecting Bi — 0, Ji = -wzri^/ii + 1)(-^ ~ 2 — i) 
reproduces the scheme used in (9| , which invokes a map- 
ping to the Jx rotation matrix of a spin ^{N — 2) particle. 
In the Toric code model, the transfer |0) to |A^ — 2) cor- 
responds to the transfer C/q \i') (the ground state with a 
single error) to Un-2 IV'} (the ground state with a logi- 
cal rotation and a single error). Thus, in a time ~ N/5 
an adversarial error model is guaranteed to propagate a 
single error into a logical gate operation and data stored 
using the Toric code does not reliably maintain its in- 
tegrity over longer time scales. 

The more widely applicable technique of Karbach and 
Stolze [lO] can be introduced to select the values of Ji 
and Bi serves to motivate the generality of the approach. 
This starts from the M x M effective Hamiltonian (M ~ 
poly(7V)), which, assuming some valid assignment for the 
Ji and Bi, can be diagonalized to find the eigenvalues Xi 
and eigenvectors |Ai). The fidelity of state transfer is 
given by 



F 



M 



-iXit 



where Ui = {M\Xi) (Ai|l). The maximum achievable 
value, for fixed {oi}, is F^ax < X)j kil- Moreover, for 
symmetric tridiagonal Hamiltonians (Ji = Jn+i-i 7^ 0), 
(Af |A,) = {-If (A,|l), so F^ax = E. I (A.|l) ? = 1. By 
selecting new eigenvalues such that, for some t and 9, 



-iXit 



'sign(aj) 



(1) 



equality can be achieved for -Flnax- To specify the eigen- 
values Ai, one selects a time t 3> tt/A' to be the state 
transfer time, where A' = Tuini^j \Xi — Xj\. Each of the 
Xi is truncated to an accuracy of 7r/i such that the num- 
ber of intervals k of 7r/t that separate A^ and Xj deter- 
mines the value of e^*^'*''^'^^^* = (—1)*^. Each eigenvalue 
can then be shifted by no more than 1 interval to ensure 
that Eqn. ([T]) is satisfied. Finally, for symmetric effective 
Hamiltonians, an inverse eigenvalue problem (lEP) |11| 
can be solved to match the new eigenvalues, and per- 
fect state transfer is assured |l5|. The linearity of the 



problem 11] guarantees that the Ji are only shifted by 
a small amount, 0(1/ ((5i)), from their initial values and 
hence the Hamiltonian perturbations are still precisely 
that - perturbations satisfying \\hi\\ < 1. For example, 
using the Toric code with Ji — ^, A' ^ 5/N'^, and the 
logical rotation occurs in a time ~ N"^ /5. 

Having designed the Hamiltonian perturbation assum- 
ing that an error occurs on a specific qubit, it is worth ob- 
serving that this is not the limit of what can be achieved. 
For example, the construction only acts on a single row 



of the lattice, and similar terms can be implemented on 
each row (the terms from one row commute with those 
on other rows). Furthermore, if an error were to occur 
on a different qubit in the row, then [8| reveals that this 
is propagated to a single error on a different spin in the 
state transfer time. This is proved by using a duality 
mapping to the perfect state transfer chain 9] , which im- 
plements t he p erfect mirroring of states in all excitation 
subspaces [l^, due to the applicability of the Jordan- 
Wigner transformation. Denoting the controUed-NOT 
gate, with control i and target j, by Cj, the mapping 
can be given as 



where 



V ■ 



^N-2 



+ 1,-1)^(21+1,1) 



(2J,0) 



'a 



(2»,0) 



Ylc^^^ 



2,0) 



The order of the second product is important, because 
the controUed-NOTs do not commute; the right-most 

term should be Wool ^^'^' ^^^^^ ^^ ^^^ ^'"^^ S^^^ ^° ^^ 
applied) and i increases to the left. The first excita- 
tion subspace of the transformed Hamiltonian maps to 
{Ui IV')}; E^nd the second excitation subspace maps to 
{UiUj |V')}j>i- Hence, any single error at a site (2«,0) 
with I > is described by UiUi-i IV') and transfers 
to UN^i-iUN-i-2 IV') in the perfect state transfer time. 
During the evolution between these two states, there is 
some probability that it has gone through an intermedi- 
ate state from which the stored information cannot be 
reliably recovered. 

Self Correction: The result on the Toric code is appli- 
cable in a wide range of cases, since the lEP technique 
allows almost any effective Hamiltonian to be recast into 
a perfect transfer system (provided that a single error 
can be connected by a polynomial sequence of local op- 
erations to a state that gets falsely corrected, as is cer- 
tainly the case for stabilizer Hamiltonians), giving a high 
likelihood of destruction of the stored data. This state- 
ment is predicated on the assumption that the minimum 
gap between eigenvalues, A', is not exponentially small. 
We shall now analyze the 2D Ising model, the archetypal 
construction of a self-correcting code for classical data, 
to see how it coincides with these conditions. Note that 
proving that the 2D Ising model is protected against this 
class of noise does not constitute a proof that the system 
is self correcting. Rather, the proof can only be used to 
show that certain systems are not self-correcting. 

The Ising Hamiltonian is defined as 



Hj 



-m 



Zi Zin 



-.3) 



where the sum is over nearest neighbour pairs, (i,j), of 
an N X N square lattice with periodic boundary con- 
ditions and qubits are placed on the vertices. The two 



ground states |0) and |1) are capable of storing 

classical information (a single Z error destroys phase in- 
formation). If a pattern of X errors affects one of the 
states, the resultant state is an eigenvector of Hj, and 
has an energy equal to the surface area of the pattern. 
To convert between the two ground states requires a pat- 
tern of X rotations on every single qubit. One can con- 
sider writing this as a progressive sequence, starting from 
a single X error, adding one rotation at each step, pro- 
gressively filling up entire rows (this sequence provides a 
convenient labelling of the qubits, 1 to N'^). All other 
choices of sequence have a similar behaviour. Since the 
stabilizers of this Hamiltonian are different to those of 
the Toric code, SH requires revision: 



JV^-2 



5H=^6J2 J^X^+l{1-ZA+2) + \5Y,B^{1-Za^+l)■ 



Motivation for the design is identical, and coincides with 
Q. To ease the analysis, without affecting the eventual 
outcome, 5H can be further modified - in the terms which 
would introduce an X which would complete an entire 
row [Xisii for integer i), we instead apply XjsaXi^i+i, 
which also starts the next row, keeping the surface area 
of the block of Xs constant. Taking Ji — \,Bi — Q, the 
effective Hamiltonian (which is an M x M matrix, where 
M = N{N - 1) - 2) can be written as 

M-l 



ifeff ={N+l)\ + 5Y, {\i) {i + l\ + \i+ 1) 

4=1 

N-1 

- ^ 2(iV - i)(\i) {i\ -t- |M -h 1 - i) (M -h 1 - 



The first step is to determine the eigenvalues. Since 5 is 
small, perturbation theory can be applied to the initial 
eigenvalues of 2(i -I- 1) for i = 1 ... TV. All of these are 
repeated twice, except for 2(A'' -|- 1), which is repeated 
(A^— 1)(A^ — 2) — 2 times. Degenerate perturbation theory 
is first applied to the repeated 2(A'' -|- 1) energy, which 
splits the degeneracy at first order. 



Ei ^ 2{N + I) + 2S cos 



{N -1){N ~2)- 1 



for integer I < z < (A^ — 1)(A^ — 2) — 2. Now, considering 
all the other energies 2(i -|- 1), by symmetry, all possible 
shifts for the repeated eigenvalues are identical, until the 
(M -|- 1 — 2i)*'' order. So, for i = 1, the energy difference 
between the two originally degenerate levels is 0{5^ ), 
and the strategy of fixing the eigenvalues would require 
an exponentially long evolution, which is not prohibitive. 
In fact, no set of quasi- local perturbations can give de- 
struction of the data; no Hamiltonian term acting locally 
on k qubits can change a state Ui jV') to any other state 
Uj IV-') for |i — j| > fc, so the effective Hamiltonian is lim- 
ited to being fc-banded, i.e. only the first k sub-diagonals 
can be non-zero. Furthermore, these are all of strength 



< 5. Thus, the perturbation calculation of eigenvalue 
splittings applies equally to all possible perturbations. 
By identical arguments, one can similarly prove that the 
4D Toric code 2*] is stable for storage of quantum infor- 
mation with respect to this noise model. 

The critical feature of the 2D Ising model is the fact 
that the energies of the sequence of X operations changes 
with the length of the sequence, scaling with N, which 
is the property known as string tension. This is essential 
because our technique has to be able to protect against 
all local errors, not just the single error that has been 
considered so far. In order to achieve this, it is nec- 
essary that when considering (degenerate) perturbation 
theory on the effective Hamiltonian, in order to connect 
any pair of initial and final states, the number of inter- 
mediate steps must grow with N and each sequence of 
steps between unperturbed eigenstates must go through 
a number of distinct energy states that grows with N, 
necessitating an increasing order of perturbation theory. 
Finally, in order for the perturbation to remain a pertur- 
bation, the smallest gap between unperturbed eigenval- 
ues must be constant. So, the spectrum for the string 
must vary with length, and the maximum energy must 
grow with N. 

In Conclusion, we have given an explicit error mech- 
anism that shows that the 2D Toric code is not self- 
correcting i.e. that stored data survives for no longer than 
0{N), where there are N^ qubits in the system [l^. This 
mechanism, in which Hamiltonian perturbations play a 
leading role, if not the precise scaling, is capable of de- 
stroying the information stored in many systems that 
do not exhibit string tension, motivating the expecta- 
tion that string tension is a necessary condition for self- 
correcting quantum memories. In contrast to previous 



approaches, we have been able to consider a dynamical 
process that does not rely on equilibrium properties, but 
describes how long storage is possible for on approach 
to equilibrium. Note, however, that our treatment has 
little to do with the thermalization process itself due to 
the limited consideration made of the bath interaction, 
as indicated by the high purity preserved in our model. 

The interplay between noise and Hamiltonian pertur- 
bations is a feature that may prove important in a num- 
ber of other settings. For example, in topological quan- 
tum computation, a non-Abelian model is defined by a 
Hamiltonian on a 2D lattice, and the idea is to store and 
compute using the excitations of the model by braiding 
them around each other. A fault occurs due to noise if a 
pair of anyons is created by some noise event and braids 
around one of the computational anyons. Unlike the case 
of memories, in topological quantum computation, we 
have the ability to regularly perform active error correc- 
tion. The general expectation is that the fault-tolerant 
threshold should be similar to that for the circuit model 
of computation, but that the inherent robustness of the 
gates should make it easier to get the gate fidelities over 
the required threshold. Since perturbations in the Hamil- 
tonian can provide the propagation mechanism to convert 
localized faults into genuine gate errors in the topological 
model, this implies that in schemes without string ten- 
sion between the anyon pairs, the expected robustness of 
the braiding operations is reduced. 
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